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The drift acceleration due to radiation reaction for a single electron in an ultraintense plane wave 
(a = eE/racw ~ 1) of arbitrary waveform and polarization is calculated and shown to be proportional 
to in the high-a limit. The cyclotron motion of an electron in a constant magnetic field and an 
ultraintense plane wave is numerically found to be quasiperiodic even in the high-a limit if the 
magnetic field is not too strong, as suggested by previous analytical work. A strong magnetic field 
causes highly chaotic electron motion and the boundary of the highly chaotic region of parameter 
space is determined numerically and shown to agree with analytical predictions. 
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It has been known for many years that qualitative 
changes occur in the behavior of an electron moving in 
an electromagnetic plane wave when the dimensionless 
strength a = eE/mcuj is of order unity. Recent ad- 
vances in laser pulse compression and amplification |^] 
have made it possible to attain such ultraintense waves 
in the laboratory and led to new investigations of their 
properties. One important effect is that an electron mov- 
ing in an ultraintense electromagnetic plane wave and 
also subjected to slowly varying "background" fields be- 
haves approximately like a particle of enhanced mass 
77170 — m^Jl + e^{A^Ai^) /mP-d^ drifting in the back- 
ground fields. Here A is the vector potential and 70 = 
\/l + a^/2 for a linearly polarized monochromatic wave. 
The fast motion in the wave can be rigorously averaged 
over if the background fields are sufficiently weak and 
the plane wave is not so strong that pair creation effects 
become significant. 

The first part of this Letter uses the guiding-center 
equations obtained in the derivation of the enhanced- 
mass approximation |^] to calculate the drift acceleration 
of an electron in a plane wave due to the radiation damp- 
ing and reaction forces. In a strong wave, the variation of 
the electron effective mass 7777 over the wave orbit must 
be properly accounted for to obtain the correct drift ac- 
celeration. The average drift acceleration is the quantity 
of greatest interest because experiments and astrophysi- 
cal phenomena typically involve many wave periods. The 
guiding-center equations provide the tools to carry out 
the averaging and determine the motion of electrons un- 
der the (previously calculated) Lorentz-Dirac radiation 
force 1^,^ . The second part of this Letter considers the 
destruction of the enhanced-mass behavior and transi- 
tion to stochasticity as the strength of the background 
fields is increased to violate the field-strength condition 
ei^back/^^cwwavc ^ 1 ucccssary for the enhanced-mass 
derivation. The breakdown of the enhanced-mass picture 
is shown numerically to predict the onset of stochasticity 
even for very strong wave intensity. Both the nonlinear 
radiation acceleration and the stochasticity are expected 



to be significant in astrophysical situations, and the first 
effect may well be visible in laboratory experiments with 
ultraintense laser pulses. 

Classical calculations are valid for high a as long as 
the discrete photon energy fiuj <^ mc^ and the amplitude 
for QED processes such as e+e^ pair creation is small 
{eEfi/mc <C 777c^) |5|. The classical high-a regime in- 
cludes a number of existing and proposed accelerator de- 
signs, such as the plasma wakefield and beat-wave accel- 
erators 1^. Lasers used in the National Ignition Facility 
and similar inertial fusion projects have a ~ 1 so strong- 
field effects are significant. Many astrophysical problems 
also involve high-a radiation sources, and in particular 
consequences of strong-wave radiation damping are dis- 
cussed in . 

The motion of an electron in a plane wave of arbitrary 
intensity is integrable both classically and within the 
Dirac equation. This occurs because a third conserved 
quantity j — Px exists in addition to the two transverse 
generalized momenta. Here 7 is the electron Lorentz fac- 
tor and px is the component of electron momentum along 
the wave axis. Including radiative effects or adding ad- 
ditional fields generally violates this conservation. The 
result of the scattering of the wave by the electron is a 
force on the electron along the wave axis, which has the 
form F — 2e^ (&) / im? for weak plane waves. In the 
high-a regime the electron radiates much more strongly 
(oc a'* rather than a? in the low-a limit) and radiates high 
harmonics of the wave frequency js]^. 

Previous work on the motion caused by radiation reac- 
tion in a strong wave has included numerical studies 
and solutions to the Lorentz-Dirac equation (in the Lan- 
dau approximation ) in special cases such as monochro- 
matic linearly and circularly polarized waves The 
guiding-center formalism gives a simple result for the 
drift velocity, the quantity of primary interest, valid for 
an arbitrarily polarized, polychromatic wave. The two 
approximations involved are the Landau approximation 
that the radiation force is calculated from the trajectory 
in the absence of the radiation force, and the guiding- 
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center approximation that the changes in the drift ve- 
locity over a single orbit are nonrelativistic and hence 
add linearly, which is well satisfied for radiation reaction 
in the classical regime. An additional advantage of the 
guiding-center formalism is that multiple applied forces 
effectively superpose, which is not in general the case for 
nonlinear equations of motion. 

The motion in a wave along x with dimensionless 
vector potential A{t — x/c) = A{ri) is given in the 
"drift frame" where the electron has zero average veloc- 
ity by Pi = —mcAi for the transverse components, and 
Px = mc(A^ — (A^))/27o, with () indicating averaging 
over r] and 70^ = 1-1- (A^) . The wave envelope is assumed 
constant in 77 so that the averages are well-defined. Vari- 
ation in the wave envelope causes a ponderomotive force, 
reviewed below. 

The radiation reaction force from scattered photons 
causes the drift velocity (the velocity of the drift frame) 
to change in time. The drift velocity component along 
the wave axis changes for an applied force F (in addition 
to the force from the wave) according to 1^ 
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Note that the instantaneous drift acceleration from an 
applied force is not necessarily in the direction of the 
applied force. The causality problems with the Lorentz- 
Dirac classical radiation force |^ do not appear in the 
Landau approximation, which is valid if the radiation 
force is a small perturbation on the motion causing the 
radiation. The instantaneous radiation force on the elec- 
tron is 
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Upon substituting (g) into (|i|), we obtain the instan- 
taneous acceleration (w.r.t. phase) in the x-direction: 
dv^/dr] = 26^70 (dA/dr/)V3mc2. Now the averaging over 
r] is simple. For a wave with intensity per frequency in- 
terval the result is 
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The acceleration is thus independent of polarization and 
the phasing between different frequences, just as in the 
nonrelativistic limit, but scales with a?' rather than a? 
for large a. 

For a linearly polarized monochromatic wave, the aver- 
age acceleration is 2e^ijj'^{a? /2 + a? + / 2 / imc^ . 
For circular polarization, the result is 2e^iJ^{a^ + 
a**) Vl + f?/3mc2. These are essentially identical to the 
considerably less transparent solutions in For circu- 
lar polarization this is just the radiated power divided 



by mc, with m = m7o. For Hnear polarization this is 
not the case, but the correct acceleration is obtained if 
we account for the fact that the electron's total radi- 
ated momentum is nonzero, unlike for circular polariza- 
tion [Q. The correct average radiation acceleration is 
thus obtained from conservation of energy and momen- 
tum if a drifting electron is assigned momentum m^foVd- 
However, there is no a priori reason for this assignment 
without the use of (|l]). 

The average radiation force can be defined as the re- 
quired constant force to balance the above acceleration: 
-Fiad = '7i7oAi.ad- The radiation force from a blackbody 
source of intensity / and temperature T is 
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Changes in the envelope of a plane wave cause a pon- 
deromotive force, which conserves ^ — Px is directed for- 
ward when the wave is rising and backward when the 
wave is falling. (There are also transverse ponderomo- 
tive forces in real beams arising from the finite spot size; 
these have a different character and are not discussed 
here.) The radiation force does not conserve j — Px and 
always acts in the forward direction. Treating an elec- 
tron in a plane wave as an enhanced-mass particle acted 
upon by ponderomotive and radiation forces gives a sim- 
ple and accurate description of single-particle behavior 
in the classical regime. 

Now we turn to consider the second problem men- 
tioned in the opening: the destruction of the enhanced- 
mass picture when strong constant electromagnetic fields 
are added to the plane wave. It will be shown that for 
one typical field configuration, i.e., one without special 
symmetries giving rise to integrability, the breakdown of 
the enhanced-mass picture is consistent with the ana- 
lytical predictions of |^ and associated with the onset 
of stochasticity over a wide range of beam intensity. An 
example of an integrable configuration is an applied mag- 
netic field parallel to the wave axis In the following 
the wave will be taken to have constant amplitude and 
the radiation force will be neglected. We study a con- 
stant magnetic field B = Bz perpendicular to a linearly 
polarized wave Ay{^) traveling in the x direction. Then 
Pz is constant and we take Pz — so that the electron mo- 
tion is confined to the xy plane. For small wave strength 
a^j = eE/mcw the motion can be analyzed perturbatively 
because the equations of motion are nearly linear |pl]] . 

When a-u, is of order unity, the equations are strongly 
nonlinear and new phenomena appear. However, the mo- 
tion is still simple for ay, > 1 as long as the applied 
magnetic field is not too strong. In the derivation of the 
equations for the motion of the guiding-center, it was 
necessary to assume — eB/mcuj = uJc/oj ^ 1, i.e., 
the electron is far from resonance. Fig. Q shows a typi- 
cal trajectory in the enhanced-mass regime: the electron 
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executes fast oscillations in the wave, while its guiding- 
center makes a slow orbit in the magnetic field. The 
relativistic nonlinearity in the equations of motion can 
destroy quasiperiodicity, but only for Wc/w ~ 1. 

With no wave present the gyrocenter of the electron 
motion (a;o, yo) ~ +py/eB, y — p^/eB) (with m = c = 
1) is constant in time. The electron still has an exactly 
constant gyrocenter for an arbitrarily strong wave: 
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Such a gyrocenter exists for any direction of the magnetic 
field Bob and any polarization of the wave: Tc = r — (r • 

b)b - b X (p + eA(0 + 1 - 7k)/eBo- 

Now consider again the particular case k || x, A || y, 
B II z. Taking pz = confines the motion to the xy 
plane so that phase space is five-dimensional (position 
{x,y), momentum {px,Py), and time t). The constants 
(^ reduce the effective dimension of phase space by two. 
Choosing particular values — Ky = of the constants 
corresponds to shifting the gyrocenters of all possible tra- 
jectories to the origin, removing two translational degrees 
of freedom. After this shift, px and py are no longer in- 
dependent coordinates but rather functions of x and y 
determined by (||). 

The existence of two constants of motion reduces the 
effective phase space in this particular case from five 
dimensions to three. Hamiltonian motion in a two- 
dimensional phase space is always integrable, so that 
three-dimensional motions such as the driven pendulum 
and the Chirikov- Taylor problem jl^ are the simplest 
that can exhibit nonintegrable behavior. The equations 
of motion after a change to the independent variable 
•q = t — x/c and introducing dimcnsionless parameters 
ab — i^c/^, ciw = eEf)/mcLO^ and a monochromatic wave 
Airf) are {oj — m — c = 1): 
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These equations are integrated numerically for various 
values of a?,, and the initial conditions x{r]a),y{rio). 

The equations of motion (^ are periodic in rj and plot- 
ting the electron's {x, y) coordinates after each period 
gives an area-preserving map of the plane to itself. Fig. |^ 
shows trajectories under this map. For low values of the 
electron initial energy, the motion is quasiperiodic and 
nearly circular, while for large values of the initial en- 
ergy the motion is chaotic, as demonstrated by numer- 
ically calculated Liapunov exponents (Fig. ^). As time 
increases the largest exponent for trajectories beginning 
on points D, E, F remains positive, indicating that neigh- 
boring initial points separate exponentially rapidly JTsf . 
The electron's initial energy affects the character of the 



motion because ab can be much larger in the rest frame 
of an electron than in the lab frame if the electron has a 
large initial velocity. To eliminate this effect the initial 
electron drift velocity is fixed at 0.5c henceforth. 

The enhanced-mass description predicts that the 
guiding-center of the electron moves in a circle: 
(xgc(t),?/gc(i)) = jxp + jvd/^) sinflt, yo + [vd /^) cosj lt), 
with Vd — 0.5c and — LOc/lold = abOJ^/l — Wd'^/To- 
At the end of each wave period, Xgc{t) and ygdt) are 
compared with the actual location of the guiding-center 
calculated from (^. As a dimcnsionless error measure 
we use the normalized sum of squares error over a cy- 
clotron orbit: Eg^ = {yt/vdf J2 (r(i) - rgc(^))^- The er- 
ror is found to be quite small {Eg^ < 0.01) for all values 
of studied as long as < 0.04. For each value of 
aw, Egc increases rapidly to order unity once ab reaches 
a certain critical field: as a threshold we define a™*(au,) 
as the value of ab where Egc = 0.01. Consistent with the 
predictions of the enhanced-mass picture, a™* remains 
nonzero for large a^ and in fact increases slightly with 
aw Even though high makes the equations (^ quite 
nonlinear, the motion remains quasiperiodic and nearly 
circular for ab < a"'*. 

As ab increases above a™*, the trajectory with initial 
velocity 0.5c becomes chaotic (has a positive Liapunov 
exponent) at some value a^. For large aw, above al the 
motion is strongly chaotic and the electron energy fluc- 
tuates wildly. This differs from nearly linear resonance 
at small a^ in that no tuning of frequencies is necessary 
for energy gain and hence energy gain is not limited by 
relativistic detuning. Fig. ^ shows numerical curves for 
a"'' and a'l as part of a schematic phase diagram. The 
two curves are adjacent over four decades of beam in- 
tensity. The dotted line (which was not calculated and 
is only schematic) separating the quasilinear resonance 
phase from the strong stochasticity phase can be defined 
by the destruction of the last invariant torus at large en- 
ergy, since in three dimensions such a torus bounds the 
energy of trajectories contained within it. Rax has previ- 
ously proposed that the stochastic motion of electrons in 
multiple plane waves may give rise to high-energy cosmic 
rays jlj] . The considerations above indicate that a single 
plane wave, together with a sufficiently strong magnetic 
field, is sufficient. 

We verified that approximately the same boundary for 
the guiding-center description applies when the wave con- 
tains two or three frequencies. It seems natural to conjec- 
ture that the guiding-center region in Fig. ^ also describes 
waves of finite bandwidth and other orientations of the 
magnetic field (excluding the integrable case B || k). 

The author wishes to thank Deepak Dhar for many 
helpful conversations. This work was supported by a U.S. 
Fulbright grant and a fellowship from the Hertz Founda- 
tion. 
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FIG. 1. Typical electron cyclotron motion with aw = 1.5, 
at = 0.02, B I z and k || x. 




FIG. 2. Surface of section showing trajectories from six 
different initial conditions with = 0.1, a;, = 0.18. 
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FIG. 3. Numerical largest Liapunov exponent ai calculated 
at different times for trajectories starting on the six labeled 
points in Fig. 2. 
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FIG. 4. Numerical values a"'* and for various Ou, as 
part of schematic phase diagram. Error bars are shown for al 

because it is difficult to determine precisely when the largest 
Liapunov exponent becomes positive. 



